Pulsar  Astrometry  by  Using  VLBI 
(PSR0329  between  Kashima  26m«5  &  Usuda  64md) 

S.  Hama,  H.Kluchi,  Y.Hanado,  Y.Takahashi,  M.lmaa 

Communications  Research  Laboratory 
€93  Hirai,  Kashima,  Ibaraki,  314  Japan 

K. Fujisawa,  H. Hi rabayash i ,  H.Kobayashi 

ISAS,  3-1  Yoshinodai,  Sagamihara,  Kanagawa,  229  Japan 

Abstract 

Communications  Research  Laboratory  (CRL),  carried  out  a  VL8 I  (Very  Lons  Baseline 
Interferometer)  experiment  of  the  pulsar  0329+54  on  November  in  1992  between 
Kashima  26md  and  Usuda  64md.  We  got  a  correlation  for  it  by  using  K-3  VLBI 
correlator,  which  was  developed  by  CRL.  Though  we  observed  a  slow  pulsar  this 
time,  we  are  going  to  make  millisecond  pulsar  experiments  by  making  good  use  of 
this  result  and  by  using  K-4  correlator  which  is  under  developing  in  CRL.  Its 
result  will  give  us  precise  positions  and  proper  motions  of  pulsars  which  are 
useful  for  obtaining  precise  pulsar  timing. 

Introduction 

Many  radio  observatories  are  interested  in  measurement  of  precise  pulsar  timing. 
Its  position  and  the  proper  motion  are  important  parameters  to  estimate  the 

pulsar  timing.  There  are  two  major  methods  to  get  them - deriving  from  the 

timing  data  itself,  and  VLBI.  The  problems  of  each  method  are  as  follows; 
timing:  precise  and  stable  data  over  a  long  range  is  necessary 
a  good  software  for  parameter  estimation  is  necessary 
VLBI:  two  large  antennas  are  necessary 

"gating"  is  necessary  to  throw  away  the  data  which  doesn't  contain 
the  pulse  itself 

Pulsar  VLBI  was  conducted  several  times  so  far.  For  example, 

by  VLA  ^  ^  -  short  baselines 

(2) 

by  VLBI  -  making  use  of  gating  function  of  Mk-llla  correlator 

by  VLBI  ^  -  gating  by  software 

The  pulsar  gating  function  of  K-3  correlator  was  used  this  time. 

Correlation  Amplitude  of  Pulssr  VLBI 

The  condition  required  for  pulsar  VLBI  is  discussed  here.  In  Japan,  as  Nobeyama 
45md  antenna  is  dedicated  for  frequency  over  10GHz,  Usuda  64md  and  Kashima 
34md  are  taken  as  the  best  couple  of  antenna.  In  practice  Kashima  34md  has 
been  in  trouble  since  this  August.  Kashima  26md  antenna,  therefore,  is  select¬ 
ed  for  this  estimation.  Let  them  have  suffix  ”u"  and  "k".  Then  estimated 
correlated  amplitude  P  is; 

eDxDuS  ./liL^xiO-26  (1) 

8  k  V  T kTu 

where  D;  antenna  diameter  [ml 
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;  antenna  efficiency 
T;  system  temperature  [K] 
s;  flux  density  of  the  star  IJy] 
k;  Boltzmann  constant  CJ/K] 

Generally  speakin9  data  of  pulsar  strength  and  profile  is  not  so  sufficient. 
One  of  the  reason  is  they  are  not  stable  over  timB  and  frequency.  So*  this  is 
just  a  rough  estimation.  In  the  case  of  PSR0329.  Downs  et.al.^*  gives  its  peak 
amplitude  in  2388MHz  as  2.6Jy  and  Lyne  et.al.^  gives  its  profile  as  shown  in 
Fig.1.  Then  we  can  adopt  an  average  strength  0.5Jy  over  50ms.  By  using  the 
following  parameters  ; 

Dk=26Im],  Du=64[m],  >7k=0.5.  Tu»0.7.  Tk=100[Kl,  TU=30[K1.  k=1.38E-23 
p  =  1. 79X  10-4  (2) 

is  obtained.  This  is  not  easy  to  detect  for  100-200  second  integration  time 
which  is  popular  for  geodetic  VLB1  . 

So  it  is  useful  to  gate  out  the  data  which  doesn't  contain  pulse  and  to  inte¬ 
grate  only  the  pulse  as  shown  in  Fig. 2.  As  it  is  generally  hard  to  know  the  a 
priori  position  of  the  pulse  in  the  period,  the  gating  positions  are  set  n 
different  wavs  using  n  correlation  units.  Gating  by  software  is  also  possible 
and  is  tried  successfully  by  Petit  et.al. 

How  gating  improves  correlated  amplitude  and  S/N  ?  Let's  think  about  the  case 
shown  in  Fig. 2  where  noise  has  average  amplitude  1  and  a  rectangle  pulse  signal 
has  amplitude  A  and  duration  1/n  of  the  pulse  period. 

In  a  normal  correlation  processing.  S  (signal)  is  A/n  and  N  (noise)  is  1,  which 
makes  S/N  =  A/n 

In  the  case  of  gating  the  pulse,  the  integration  time  is  1/n  which  makes  the 
average  amplitude  of  the  noise  *Jr\.  As  S  is  A,  S/N  becomes  AA/fT.  Thus  S/N  is 
improved  -/rT  times.  In  case  n=8  and  the  antenna  pair  is  the  Kashima  26md  and 
the  Usuda  64md.  the  detection  of  the  fringe  is  possible  as  follows; 

P  =  l. 43X10“3  (3> 

It  is  true  that  n  times  long  observation  also  improves  its  S/N  Jn  times,  but 
gating  has  the  following  merits: 

1.  Restrictions  of  the  antenna  machine  time  and  slewing  are  less 

2.  Not  sensitive  to  the  long-term  instability  of  the  frequency  standard 

3.  Tape  consumption  and  restrictions  of  the  tape  length  are  less 

4.  Good  for  faint  pulsars  because  the  correlated  amplitude  itself  is  l?rge 

Correlation  Processing  System  of  CRL  and  Its  problems 

The  correlation  processing  system  of  CRL  is  shown  in  Fig. 3.  Both  K-3  and  K-4 
tape  are  possible  to  be  processed.  The  software  NKROSS  was  modified  for  pulsar 
processing.  It  is  written  with  HP  Basic  and  run  on  a  personal  computer  HP330 
(the  CPU  is  68020).  A  priori  parameter  is  given  from  the  HP330  to  the  correla¬ 
tor  and  the  correlation  data  is  acquired  by  the  host  through  GP I B  every  PP 
(parameter  period;  usually  one  through  four  seconds).  The  gating  function  of 
the  K-3  correlation  processor  is  only  once  per  PP.  So  we  selected  PSR0329+54 
which  has  long  period  as  714.5ms  and  is  strong  enough  to  get  the  fringe.  The 
ideal  PP  should  be  the  same  as  the  pulse  period.  But.  as  PP  is  quantized  by  5ms 
we  set  PP  1  second  which  is  close  to  the  pulse  period.  The  timing  of  opening 
and  closing  of  the  gate  can  be  set  in  the  unit  of  bit  (250ns)  by  the  host. 
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The  dispersion  is  always  problem  for  pulsar  observation.  The  pulse  period  is 
stretched  out  to  At  as^'; 

At=  8.3E3*DM*RF~3*BW  tsecl  (4) 

where,  DM:  Dispersion  Measure  [pc/cm3] 

RF :  Observing  Frequency  [MHz] 

BW:  Bandwidth  [MHz] 

In  this  experiment  PSR0329  gives 

DM=27 [pc/cm3] ,  RF=2300 [MHz] ,  BW=2[MHz]  then, 

A  t=  36. 8//  s  (5) 

is  obtained.  This  value  is  far  smaller  than  50tns  which  is  estimated  as  the 
pulse  duration  in  the  last  section.  Moreover,  when  the  bandwidth  ranges  over 
100MHz  for  bandwidth  synthesis.  At  is  only  1.97ms  which  is  still  smaller  than 
the  50ms.  Of  course  this  should  be  compensated  in  the  case  of  bandwidth  synthe¬ 
sis. 

A  small  problem  appeared  in  the  actual  processing.  Fringe  rotation  (=RF*<^T, 
dX  is  delay  rate)  was  very  small  because  the  baseline  is  short  (200km),  frequen¬ 
cy  is  low  (2.3GHz),  the  right  ascension  is  high.  Our  K-3  correlator  makes  use 
of  3  I  eve  I -qu ant i zed  sinusoidal  function  (as  Fig. 4)  to  compensate  fringe  rota¬ 
tion.  On  the  other  hand  the  data  to  be  integrated  within  a  PP  is  small  because 
of  gating.  Some  PP’s,  therefore,  have  most  data  on  the  co r re  1  a t i on- r e s t r i c t ed 
area.  It  means  that  such  PP  is  given  heavy  weighting  even  if  its  quality  is 
bad.  This  is  notable  in  the  case  of  short  period  pulsars. 

The  Result  of  the  Experiment 

The  parameters  of  the  experiments  are  shown  in  Table  1.  Strong  radio  source 
3C84  is  used  for  an  initial  clock  search.  The  period  P=714.5ms  is  divided  by  n 

so  that  each  gating  time  is  1 /n  of  the  pulse  period. 

Fig. 5  shows  the  result  when  the  data  is  processed  in  normal  geodetic  mode. 
Fig. 6  is  processed  with  the  gating  of  P/16  over  the  same  integration  time. 
Fig. 6  has  a  clearer  fringe  than  Fig. 5  because  the  S/N  is  improved  as  mentioned 
before.  Sometimes  PP  with  a  large  correlation  amplitude  is  seen  in  Fig. 6.  The 
reason  is  probably  that  mentioned  in  the  last  section. 

Table  2  shows  the  correlation  amplitude  for  both  gated  to  period/16  and  nan- 
gated  when  the  parameter  is  integration  time.  As  there  is  only  one  gating  (P/16 
in  this  case)  in  a  PP  (=  1  sec  this  time),  the  efficiency  of  the  integration  time 
is  P/1.  Gating,  therefore,  improved  the  amplitude  by  the  factor  1 6*P/1 (=1 1 .4) . 

In  Table  3  the  amplitude  is  compared  between  pulse  area  and  non-pulse  area. 

In  Table  4  the  parameter  is  n;  the  number  the  period  is  divided  with.  The 
integration  time  is  480sec.  The  amplitude  increases  almost  in  proportion  to  n. 

Pulsar  VLBI  for  Millisecond  Pulsars 

We  plan  to  make  VLBI  experiment  also  for  millisecond  pulsar.  In  the  case  of 
1937+21,  let's  confirm  if  we  get  its  fringe  by  using  E  q .  1  .  Kashima  34md  is 
taken  this  time  instead  of  the  26md  antenna.  As  shown  in  Table  5,  average 
intensity*65  is  so  small  that  normal  correlation  gives  the  fringe  which  is  under 
the  limit  of  the  detection.  But  after  the  gating  of  P/16(=97us),  fringe  can  be 
detected.  The  dispersion  of  1937+21,  whose  DM  is  71 [pc/cm3],  becomes; 
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At  =  43.7ns  for  RF=1.5GHz,  BW=250kHz 

At  =  349.2ns  for  RF=1 .5GHz,  BW=2MHz 

At  =  96.9ns  for  RF=2.2GHz,  BW=2MHz 

In  the  case  of  1.5GHz,  wide  bandwidth  as  2MHz  makes  wave  wider  than  P/16.  So 

wider  gating  and  narrower  bandwidth  should  be  selected  for  lower  frequency. 
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20  40  60<ms)  Normal  S  =  A/n  N  =  1  S/N=A/n 

Gated  S  =  A  N  =  Vn  S/N»AA/n 

Fia.1  Profile  of  PSR0329+54  Fig. 2  Gating  for  Pulsar 

Correlation  &  Integration 
only  within  the  shadowed  area 


for  ANALYSIS 


Fig. 3 


Correlation  Processing  System  of  CRL 


Table  1  Parameter  of  the  pulsar  VLBI  experiment 


Basel  i n e 

Kash ima26md-Usuda64md 

Observed  Star 

PSR0329+54/  3C84 

Frequency 

2222 ,99MHz 

Bandw i d  t  h 

2MHz 

Polarization 

RHCP 

Correlation 

K-3  correlator  at  Kashima 

correlation  restricted  area 


F  i  g  .  4 

3-level  sin  function  for  compensation  of  fringe  rotation 


Fig. 5  Normal  correlation  of  PSR0329 

Cat  2.3GHz,  BW=2MHz/  I n t e g r a t i o n=480s e c ) 


Fig. 6  Gated  correlation  of  PSR0329 

(at  2.3GHz/  BW=2MHzx  I nteg r at i on=480sec) 
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Table  2  Correlation  Amplitude 
(Gate  time  =  period/16) 


Gated 

Non-gated 

Integration  time 

4.44E-3 

0 . 35E-3 

490sec 

4.77E-3 

0.39E-3 

250sec 

4.47E-3 

0 . 46E-3 

1 25sec 

Table  3  Correlation  Amplitude 
(Integration  time  =  250sec) 


1 . 04E-3  (area  A) 
P  4.77E-3  (area  B) 
16  1  . 22E-3  (area  C) 


ABCDEFGH  IJKLMNOP 


Table  4  Correlation  Amplitude 
(Integration  time  =  480sec) 


0.35E-3  without  gating 
0.93E-3  ( n  =  4) 

2.01E-3  (n=8) 

4.44E-3  ( n= I  6 ) 


Table  5  Correlation  Amp.  of  millisecond  pulsar  1937+21 
(assumed  gating  covers  3/4  of  the  total  flux) 


Average 

1  ntens  i  ty 

W  i  t  hou  t 
gating 

Gated 
<n  =  8) 

Gated 
(n= 1 6) 

1 . 5GHz 
2.2GHz 

9m  J  y 

3 . 3m  J  y 

9 . 5E-5 
2.9E-5 

5 .7E-4 

1  .8E-4 

1  .  IE-3 

3 .5E-4 

(gating  time  i 

s  per i od /n) 

Dispersion  At=  43  7ns  for  rf=1.5GHz,  BW=250kHz 
A t  =  349.2ns  for  RF=1.5GHz,  BW  =  2MHz 
At=  96.9ns  for  RF=2.2GHz,  BW=2MHz 
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QUESTIONS  AND  ANSWERS 

D.  Allan,  Allan  Time:  Do  you  have  any  stability  measurements  against  any  reference  clocks? 
eg:  The  arrival  time  stability  against  a  reference  clock. 

S.  Hama,  CRL:  Now  we  use  a  reference  clock;  a  cesium.  The  Cesium  is  compared  to  UTC 
Japan  by  using  GPS. 

D.  Allan:  Do  you  have  some  results? 

S.  Hama:  For  timing  measurement,  I  do  not  have  the  figure,  but  we  are  conducting  timing 
measurement  experiment.  We  have  some  results,  so  I  can  show  you  later. 

D.  Allan:  Thank  you. 
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